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GLOBAL ESTIMATES FOR GREEN'S MATRIX OF SECOND ORDER 

PARABOLIC SYSTEMS WITH APPLICATION TO ELLIPTIC SYSTEMS IN 

TWO DIMENSIONAL DOMAINS 

Q ■ SUNGWON CHO, HONGJIE DONG, AND SEICK KIM 

Abstract. We establish global Gaussian estimates for the Green's matrix of divergence 
form, second order parabolic systems in a cylindrical domain under the assumption that 
weak solutions of the system vanishing on a portion of the boundary satisfy a certain local 
C 3 boundedness estimate and a local Holder estimate. From these estimates, we also derive 

f^ , global estimates for the Green's matrix for elliptic systems with bounded measurable co- 

efficients in two dimensional domains. We present a unified approach valid for both the 
scalar and vectorial cases and discuss several applications of our result. 

Ph 

C^ , 1. Introduction 

Fundamental solutions of parabolic equations in divergence form with bounded mea- 
surable coefficients have been a subject of research for many years. The first significant 
step in this direction was made in 1958 by Nash 1341 . who established certain estimates of 
K^ ' the fundamental solutions in proving local Holder continuity of weak solutions. In 1967, 

O^ i Aronson fTl proved two-sided Gaussian estimates for the fundamental solutions by using 

the parabolic Harnack inequality of Moser ll33l . In 1986, Fabes and Stroock 1151 showed 
that the idea of Nash could be used to establish Aronson's Gaussian bounds, which conse- 
quently gave a new proof of Moser's parabohc Harnack inequality. There are many books 
and articles related to this subject; see e.g., ifTOl [35l [36l [38l and references therein. 
(^ ' Compared to a long history in the study of fundamental solutions for parabohc equations 

with real coeflicients, there has been relatively little study on the fundamental matrices 
for parabolic systems until recently, except when the coefficients are sufficiently regular; 
see e.g., Eidel'man lfT4l . In 1996, Auscher ||2l gave a new proof of Aronson's Gaussian 
^\ ' upper bound for the fundamental solution of parabolic equations with time independent 

^ . coefficients, which carries over to the case of L°°-complex perturbations of real coefficients. 

We note that a parabolic equation with complex coefficients is, in fact, a special case of a 
system of parabolic equations. Since then, there has been active research in this direction; 
see e.g., i3i 4, 5, 6, 8, 22 , 25il for related results. 

In this article, we study the Green's matrix for parabolic systems 

in n 

(1.1) du'ldt - X ^ £*ff (^ 'f (f. x)Dp uJ), D„ = d/dx„, /=!,..., m, 

in a cylindrical domain 1/ = R x Q, where Q is a (possibly unbounded) domain in R". 
We assume that the coefficients are measurable functions defined in the whole space R"^' 
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satisfying the strong parabolicity condition 

m n m n 

and also the uniform boundedness condition 

m n 

(1.3) YuYMt^^^A -"'^^ S{t,x)eW^\ 

for some constant v 6 (0, 1]. We point out that the coefficients are assumed to be neither 
time independent nor symmetric. We will later impose some further assumptions on the 
operator but not explicitly on its coefficients. 

We are interested in the following global Gaussian estimate for the Green's matrix 
0{t, X, s, y) of the parabolic system dl.lb in a cylindrical domain 1/ = R x Q; For any 
r > 0, there exists a positive constant C such that for all f, s G R satisfying < t - s < T 
and X, y € Q, we have 

(1-4) mt,x,s,y)\ < '^ expj-^^^::^ 

(t - s)"i- y t - s 

where /c is a positive constant independent of T . We are also interested in the following 
global estimate for the Green's matrix 0{t, x, s, y) of the parabolic system ( 11.11 ) in 1/ = 
R X Q when the base Q. is subject to a certain condition: For any T > Q, there exists a 
positive constant C such that for all t,s eR satisfying < t - s < T and x,y € Q., we have 

(1.5) \g{t,x,s,y)\<cllA , ^^" ) (lA 



Vf- s V\x-y\) \ Vf- s V \x - y\ , 

.(t-s)-"'^cxp{-K\x-y\-/(t-s)}, 

where k > and ju G (0, 1] are constants independent of T, and we used the notation 
a /\b - min(fl, b),a\/ b - max(a, b), and d^ = dist(x, dQ). 

The goal of this article is to present how one can derive the global estimate like ( |1.4t 
for the Green's matrix using a local boundedness estimate for the weak solutions of the 
parabolic system vanishing on a portion of the boundary; see Condition (LB) below for 
the precise statement of the local boundedness estimate. In fact, we show that such a local 
boundedness estimate is a necessary and sufficient condition for the Green's matrix of the 
system to have a global estimate like jlAj . We shall also show how to derive a global es- 
timate like (11.5b for the Green's matrix by using a local Holder continuity estimate for the 
solution of the system vanishing on a portion of the boundary; see Condition (LH) below 
for the statement of the local Holder estimate. There is a standard method in constructing 
the Green's matrix for elliptic systems in a domain Q. out of the "Dirichlet heat kernel" of 
the elliptic system, namely by integrating it with respect to f-variable from zero to infinity. 
By utilizing the above global estimates ( |1.4b and ( |1.5b . we obtain the following global esti- 
mate for Green's matrix G{x, y) for elliptic systems with bounded measurable coefficients 
in a two dimensional domain fl: 

(1.6) |G(x,3;)| <c(l A-—^) (l A— ^1 |l+ln+' ^ 



\x-y\l \ \x-y\l [ \\x-y\) 

where C > and fi e (0, 1] are constants depending on Q., and we used the notation 
In+f - max(lnr, 0). We do not consider Green's matrix for elliptic systems in a three 
or higher dimensional domain in this article. For treatment of such cases, we refer to a 
very recent article ll24l . where conditions similar to ours were introduced. In fact, our 
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conditions (LB) and (LH) are motivated by the corresponding elliptic conditions appeared 
in ll24l . We point out that most of the results in ll24l can be also obtained by following the 
above mentioned "Dirichlet heat kernel" approach and utilizing the global estimates like 
( ll.4l i and ( 11.5b . albeit it would be far more complicated to do so. The axiomatic approach 
adopted in this article is also in the spirit of 1221123)1 . 

The novelty of our work is in presenting a unifying method that establishes the global 
estimates (11.4b and ( 11.5b for the Green's function for the uniformly parabolic equations 
with bounded measurable coefficients as well as for the Green's matrix of the parabolic 
systems (II. lb . for instance, with coefficients uniformly continuous or VMO in x-variables 
and bounded measurable in f-variable, in a cylindrical domain R x O, where the base O 
is a C' domain or a Lipschitz domain with a sufficiently small Lipschitz constant. Also, 
this article provides a unified approach in establishing the global estimate (11.6b for both 
Green's function for elliptic equations and Green's matrix for elliptic systems with bounded 
measurable coefficients in two dimensional domains. In a recent article lU, we proved 
the existence of Green's matrix for the parabolic system ( II. lb in an arbitrary cylindrical 
domain 'ZY = R x fi under the assumption that weak solutions of the systems satisfy an 
interior Holder continuity estimate. We also derived various local estimates for the Green's 
matrix under the same assumption but, however, the global Gaussian estimate like ( 11.4b 
was proved only in the case when Q = R". In another recent article ||121 . the second and 
third named authors proved the existence and local estimates of Green's matrix for elliptic 
systems with bounded measurable coefficients in a two dimensional domain Q with finite 
area or width, by utilizing the estimates established in |8|. However, again, the global 
estimate similar to ( 11.6b was established only when Q is the open set above a Lipschitz 
graph (^ : R ^ R. In this sense, the present article may be considered as a sequel of both 
articles ||8l and lfT2l with a considerable improvement. 

The organization of the paper is as follows. In Section |2] we introduce some notation 
and definitions including our definition of the Green's matrix of the parabolic system ( II. lb 
in 'ZY = R X fi. In Section|3] we give precise statement of the conditions (LB) and (LH) and 
state our main theorems. In Section IH we present some applications of our main results. 
The proofs of our main results are given in Section |5] and several technical lemmas are 
proved in Appendix. Section|6]is devoted to the discussion of global estimates for Green's 
matrices for elliptic systems with bounded measurable coefficients in two dimensional do- 
mains and Section |7] is allocated to a brief discussion regarding the global estimates for 
systems with Holder or Dini continuous coefficients. 

Finally, a few remarks are in order. Green's functions for uniformly elliptic equations 
wiffi bounded measurable coefficients were extensively studied in classical papers 1211 
|3T1 . Green's matrices for elliptic systems were investigated earlier in ifTSl [T6l [iTl . In 15], 
Auscher and Tchamitchian introduced the "Dirichlet property (D)" in connection with the 
Gaussian estimates for the heat kernel of the elliptic operator with complex coefficients, 
which is related to the condition (LH) of this article; see Remark [J!9] 

2. Notation and Definitions 

We mainly follow the notation and definitions of [[S], most of which were in turn chosen 
to be compatible with those used in ll28l . 

2.1. Basic notation. Let ^ be a parabolic operator acting on column vector valued func- 
tions H = (m', . . . , ul^Y defined on a domain in R"^' in the following way: 

(2.1) ^u^u,- DaiA"'^ Dpu), 
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where we use the usual summation convention over repeated indices a,/3 = !,...,«, and 
A"^ - A"^{t,x) are mx m matrix valued functions on R"^' with entries A"^ that satisfy 
the conditions ( 11.21 ) and ( 11.3b . Notice that the /-th component of the column vector ^u is 
presented in (II. 11 1. The adjoint operator 'Jif is defined by 

'^U = -U, - Dai'A"f^D/}U), 

where 'A"^ = (A^")^; i.e., A'-'j' = A^f . Notice that the coefficients A"f satisfy the condi- 
tions (11.2b and (11.3b with the same constant v. 

We use X = {t, x) to denote a point in R"^'; x = (xi,. . ., x„) will always be a point in 
R". We also write Y = (s,y), Xq = (fo,xo), etc. We define the parabolic distance between 
the points X - (t, x) and Y - {s,y) in R"+' as 



\X-Y\p:^m<ix(^J\f^\,\x-y\), 

where | ■ | denotes the usual Euclidean norm. For an open set 1/ c R"^', we denote 

dx = dist(X, dpV) = inf{|X - Y\p : Y e dpV]; inf = oo, 

where dptl denotes the usual parabolic boundary of TY. 

For a given function u - ii{X) - u(t,x) defined on 2 c R"^', we use DjU for du/dxi 
while we use u, (or sometimes DfU) for du/dt. We also write Du (or sometimes D^u) for 
the vector {D[U, . . .,D„u). For// € (0, 1], we define 

\u\^l/2,^,■,Q = |m|o;g + W^l/2,^,■,Q := sup|M(X)| + sup — — — —r^ — . 

XeQ X.YeQ \X - Y\p 

Xi=Y 



By 'T^^'^^'f'iQ) we denote the set of all bounded measurable functions m on 2 for which 
\u\fi/2,fi;Q is finite. We use the following notation for basic cylinders in 



jn+l. 



Q;(X)^(t-r\t)xBAxy, 
Qt(X)^(t,t + r^)xBAxy, 
Q,.{X)^{t-r\t + r^)xB,.{x), 

where -Sr(x) is the usual Euclidean ball of radius r centered at jc e R". For an open set 
'WcR"+', we define 

nJriX) = 1/ n Q,.{xx Sr(X) = 5pl/ n Q,.{xx 

and similarly lAfiX) and S^{X). We write tlito) for the set of all points (fo, x) in tl and 
I(!U) for the set of all t such that 'Z/(0 is nonempty. We denote 

IIImIII^ = WDuWl-,^ + ess sup \\u{t, Olliif^m,- 

re/CW) 

Finally, we denote a A b - min(fl, b) and ay b - max(a, b). 

2.2. Function spaces. Throughout this section, we shall always denote by Q the cylinder 
(a, b) X Q., where -co < a < b < oo and O is an open connected set in R". We denote by 
W^ ' (Q) the Hilbert space with the inner product 



<M,v)^o.i,g, := I Mv+V I DkiiDkV 
' Jq izi Jo 



e PT 



GREEN'S MATRIX FOR PARABOLIC SYSTEMS 



and by /C'''(2) the Hilbert space with the inner product 



'^,-\Q) 



Jq k~i-JQ -JQ 



We define 'fiiQ) as the Banach space consisting of all elements of #^ " (Q) having a finite 
norm ||m||-:^2(0 •- lll"llle- The space 1^ ''(0 is obtained by completing the set #^'"'(2) in 
the norm of fiiQ)- 

Let S c Q and m be a #^ '(2) function. We say that u vanishes (or write u - 0) onS 
if M is a Hmit in W°'\Q) of a sequence of functions in C^(Q \ S). We define /aCg) (resp. 
%°'\Q)) the set of all functions u in >2(0 (resp. f^'^iQ)) that vanishes on the lateral 
boundary d.cQ :- {a, b) x 5Q of Q. By a well known Sobolev-like embedding theorem (see 
e.g., ESI §11.3]), we have 

(2.2) l|M|lL2.4/„(e) < C(«) lllMllle, Vm € f.iQ)- 

The space ;#^ ' {tl) (I < q < oo) denotes the Banach space consisting of functions 
u G U{11) with weak derivatives D„m € WilJ) {a - !,...,«) with the norm 

We write u e U^(%1) if u e L°°(1/) has a support in 'Z7. 

2.3. Weak Solutions. For/, g^,, e Lj^^CU)"' (a = !,...,«), we say that « is a weak 
solution of JCu - f + Dag,, in 1/ if « 6 tiillT and satisfies 

Jtl J'U J%t J'U 

Similarly, we say that « is a weak solution of 'Jffu - f + Dag„ in 1/ if « e I'^CZ/)'" and 
satisfies 

(2.3) r M>i + r 'A'^fD^uJDj = r />' - r ^^A.^', V.^ e CCZ/)". 

2.4. Green's matrix. Let 1/ = R x O be a cylindrical domain, where Q. is an open con- 
nected set in R". We say that an m x m matrix valued function 0{X, Y) = 0{t, x, s,y), with 
entries gijiX, Y) defined on the set {{X, Y) ell xll : X ^ Y},is a Green's matrix of if in 
1/ if it satisfies the following properties: 

i) g(-, Y) € #:°;' {^() and ^g{-, Y) = 6yI for all 7 € 1/, in the sense that 



Jii 



2 



[lA \ QriY)) for all y e 1/ and r > and §(-, Y) vanishes on dpOA. 
iii) For any / = (/',..., /™)^ 6 L^ill), the function m given by 

u{X):^ f g(YX)f(Y)dY 
Jn 

belongs to Y^C^) and satisfies '^u - f in the sense of ( |2.3l l. 
We note that part iii) of the above definition gives the uniqueness of a Green's matrix; see 
IS). We shall hereafter say that g(X, Y) is "the" Green's matrix of Jf in 1/ = R x Q if it 
satisfies all the above properties. 
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3. Main results 



The following condition (LB) shall be used to obtain the global Gaussian estimates for 
the Green's matrix g{X, Y) = g(t, x,s,y) of ^ m^^ ^Rx Q. See TheoremlTTlbelow. 

Condition (LB). There exist /?,„„ v € (0, oo] and No > so that for sWX ell and Q < R < 
Rmax, the following holds, 
i) If H is a weak solution of J^u = in 1A^(X) vanishing on S^{X), then we have 

ii) If M is a weak solution of '^u = in t(^{X) vanishing on S^(X), then we have 

The following condition (IH) means that weak solutions of J/fu = and '^u = 
enjoy interior Holder continuity estimates with exponent jUq- It is not hard to see that this 
condition is equivalent to saying that the operator ^ and its adjoint lif satisfy the property 
(PH) in [si; see Lemma [8721 for the proof. 

Condition (IH). There exist yUo € (0, 1], /?c e (0, oo], and Co > so that for all X e 1/ and 
< R < Re /\ dx, the following holds, 
i) If H is a weak solution of Jffii = in Q]^{X), then we have 

1/2 



[uU2,,:Q,,iX) < CoR-"" 



iRpJ 






ii) If M is a weak solution of 'Jffu = in Q^(X), then we have 

( r V^' 

\JqHX) ) 

Theorem 3.1. Let tl — M^xD. and assume the conditions (IH) and (LB). Then the Green 's 
matrix QiX, Y) of ££ in 11 exists and for all X,Y etl with X ^ Y, we have 

(3.2) \@{t, X, s,y)\ < Cx(Q,.o){t - s) ■ {{t - s) A /?L.} "^' exp {-k\x - y\^/it - s)} , 

where C — C(n, m, v, Nq) and k — k(v). 

The following theorem says that the converse of Theorem B.ll is also true. 

Theorem 3.3. Assume the condition (IH). Let 0(X, Y) be the Green 's matrix of S£ in 
t( — Wx D.. Suppose there exist Rmax £ (0, o°] and positive constants Cq and k such that 
for all X,Y elA with X+Y,we have 

(3.4) \Q{t, X, s,y)\ < Co;r(0,oo)(f - s) ■ [{t - s) A RI„,\"'' exp [-k\x - y\^l{t - s)] . 

Then the condition (LB) is satisfied with the same R„,ax and Nq — Noin, m, v, Co, k). 

The following condition (LH) means, loosely speaking, that weak solutions of Jifu = 
and '^u = vanishing on a relatively open subset S of dp1( are locally Holder continuous 
up to S with exponent jUo. 

Condition (LH). There exist juo e (0, 1], R,„ax e (0, oo], and A^i > so that for all X € '1/ 
and < R < R,„ax, the following holds. 
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i) If H is a weak solution of ^u = in ^/^(X) vanishing on S^Qi), then we have 

{ r V^' 

mf,/2.^^■,Q,.,(X)<NlR-''«\-jr |h|2 ; u ^ X'U,(X)U. 
\Jqj,{X) ) 

ii) If H is a weak solution of '^u = in ^^{X) vanishing on S^{X), then we have 

( r V^' 

\Jq*{X) ) 

Remark 3.5. In the above condition (LH), the constant /?„„„ is interchangeable with cR„ 
for any fixed c e (0, oo), possibly at the cost of increasing the constant A^i . 



It is not hard to see that (LH) implies (LB); see Lemma ISTI in Appendix. Also, it is 
obvious that (LH) implies (IH). Therefore if 1/ = R x Q and (LH) is satisfied, then by 
Theorem l3.1l the Green's matrix @{X, Y) of ^ m1( exists and satisfies the estimate ( 13.21 ). 
The following theorem says that in fact, in such a case, a better bound for the Green's 
matrix 0{X, Y) is available near the boundary 5^1/ = R x dH. 

Theorem 3.6. Let 1/ = R x Q and assume the condition (LH). Let Q{X, Y) be the Green 's 
matrix of S£ in 11 and denote 



(3.7) 5(X,F)=lA lA 

I R,„a.^\X-Y\p}\ R,na.^\X-Y\, 

Then for all X,Y eV with X + Y,we have 

(3.8) \Q{t, X, s,y)\ < Cx(0.oc){t - s) ■ 6{X, Yf" {(t - s) A rL,]'"'' exp Lj.^^-^\ , 

where C — C{n, m, v,/jo, Ni) and k — k{v). 

Remark 3.9. Suppose the operator ^ satisfies the following property, which we shall refer 
to as the condition (LH'): There exist fiQ e (0, 1], R,„ax e (0, oo], and Co > such that for 
all X e 1/ and < R < Rmax, the following holds; 

i) If M is a weak solution of .if « = in tl]^{X) vanishing on >S^(X), then we have 

r |£)h|- < Co (-]" ''" r \Du\-, W<p<r<R. 

ii) If H is a weak solution of '^u = in 'ZY^(X) vanishing on S^{X), then we have 

r \Duf- <Q)(-T''° f |£>M|^ W<p<r<R. 
J'U*(X) ^f' J'uux) 

Then, the condition (LH) is satisfied with A^i = Ni{n,m,v,fj.o,Co) and the same /uo and 
R„,„/, see Lemma [84] in Appendix. The condition (LH') is reminiscent of the "Dirichlet 
property (D)", which Auscher and Tchamitchian introduced in Q in connection with the 
Gaussian estimates for the heat kernel of elliptic operators with complex coefficients. We 
note that the condition (LH) is weaker than condition (LH') in general. 



Remark 3.10. It is clear that the estimate ( 11.4b in the introduction follow from Theorem l3.1l 
Also, we note that the estimate ( 11.51 ) in the introduction follows from Theorem l3.6l if Q is 
bounded or R„,a^ - c»; see Section|4]and Section|2]for further discussion. 
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4. Some Applications of Main Results 



4.1. Scalar case. Let Q. be an arbitrary open connected set in R" and 'W = R x Q. In 
the scalar case (i.e., m - 1), both conditions (LB) and (IH) are satisfied with R„tax - 
oo and A^o = No{n,v); see e.g., 1281 [30l . Also, in the scalar case, the Green's matrix 
becomes a nonnegative scalar function. Therefore, the following corollary is an immediate 
consequence of Theorem l3.1l 

Corollary 4.1. Let 1( - R x n. Ifm = 1, then the Green' s function @{X, Y) of J£ in 1/ 
exists and for X,Y etl with X ^ Y, we have 

(4.2) g(t, X, s, y) < Cx(0.oo)(t ' ^) ■ (t - sT"'^ exp [-k\x - y\-l(t - s)] , 

where C — C{n, v) and k — k(v) are universal constants independent of Q. 

Remark 4.3. Corollary 14. H is widely known and originally due to Aronson |[T|. 

Moreover, in the scalar case, the condition (LH) is satisfied if the base Q. satisfies 
the condition (S), the definition of which is given below. In fact, if ^ is a small L°°- 
perturbation of a diagonal system, then the condition (LH) is satisfied whenever the base 
Q. satisfies the condition (S); see ^4.2l and Lemma lS^ below. 

Condition (S). There exist 6 > and Ra e (0, oo] such that 

\Br{x) \ Q| > 0\Br(x)\, Vx € dQ, V/J < R,,. 

The following corollary is then an easy consequence of Theorem l3.6l and Corollarv l4.1l 

Corollary 4.4. Let li - Rx Q., where Q satisfies the condition (S). Ifm = 1, then the 
Green 's function Q(X, Y) of J£ in 11 exists and satisfies the estimate ( 14.2b . Denote 

6{X, r) = I 1 A P- -I ( 1 A 



Rah\X-Y\p]\ RaA\X-Y\p 

Then for X,Y e 1( with X ^ Y, we also have 

g(t, X, s, y) < Cx(0,o.)(t - s) ■ 6(X, Yf" [(t - s) A R^f"'^ exp i-J^^-^\ , 

where C — C{n, v, 6), hq — yUo(n, v, 6), and k — k{v). 

Example 4.5. R" satisfies the condition (S) with 6=1/2 and Ra = oo. Then we have 

and by Corollary |431 for all X,Y ell with X ?t F, we have 

G{t, X, s, y) < Cx(Q,^){t - s) ■ 6{X, YY\t - sY"'^ exp \-k'' ~ ^' 



t- s 
where C = C(«, v), yuo = yUo(«, v), and k = k(v). 

Remark 4.6. In Corollarv l4.4l one can allow for lower order terms in ^; c.f. Corollarv l4.14l 
below. 
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4.2. L°° -perturbation of diagonal systems. Let a"^(X) be scalar functions satisfying 

n 

(4.7) a"^m^/3^. > vol^p, V^eR"; J]\a"^xf<yf, 

for all X 6 R"+' with some constant vo € (0, 1]. Let 1/ = R x Q with the base fi c R" 
satisfying the condition (S). Let A "^ be the coefficients of the operator J/f. We denote 

/ A 1/2 

(4.8) <r = sup V V \A';f{X) - a"^X)d,j\ 

where S/j is the usual Kronecker delta symbol. 

By Lemma [8^ there exists S'q - S'o{n, vq) such that if S" < <§q, then the condition (LH) 
is satisfied with jiq - jJo{n,vo,6), R„,„^ = Rg, and A^i = N\{n,m,VQ,0). Therefore, the 
following corollary is another easy consequence of Theorem l3.6l 

Corollary 4.9. Assume that a"^(X) satisfy the condition (14.7b . Let 1i -M.xQ., where fi 
satisfies the condition (S), and define 6(X, Y) as in (I3.7l l with Rmax — Ra- Let S be defined as 
in (I4.8l l, where A" (X) are the coefficients of the operator ^. There exists Sq — S'o{n, vq, 6) 
such that if S" < Sq, then the Green 's matrix Q{X, Y) of I£ in 11 exists and for all X,Y etl 
with X + Y,we have 

\g(t, X, s,y)\ < Cx(0,co)it - s) ■ 6(X, Yf" [it - s) A RI}'"'^ exp Lj.^^-^\ , 

where C, jiq, and k are constants depending on n, m, vq, and 6. 

Example 4.10. Let O = [x e R" ; x„ > ip{x')}, where x = {x',x„) and ip : R""' ^ R is 
a Lipschitz function with a Lipschitz constant K. Then Q. satisfies the condition (S) with 
- 0{n, K) and Rg - oo and we have 

6(X,Y)^\\A — 1(1 A ^ — I in'i/^RxQ. 

' ' ' \ \X-Y\p}\ \X-Y\p} 

If ^ is a small L°° -perturbation of a diagonal system in the sense of Corollary 14.91 then 
the Green's matrix 0{t, x, s, y) of ^ in 'ZY exists, and for all X,Y e^ with X ^ F, we have 



\g{t,x, s,y)\ < Cxio,o.)(t - s) ■ 6{X, YfH - i)-"/-expl-^ 



X- 



■y\' 



t - s 
where C,/io, and k are constants depending on n, m, vq, and K. 

4.3. Systems with VMO., coefficients. For a measurable function / = f(X) - f(t,x) 
defined on R"^', we set forp > 



Upif):^ sup sup -f -f \f(y,s)- f,As)\dyds\ f^,{s)^-V f(s 

XeR"+' r<p Jt-r- JBAx) J BAx) 



We say that / belongs to VMO.v if limp^o (JJpif) - 0. Note that VMO^ is a stricdy larger 
class than the classical VMO space. In particular, VMOx contains all functions uniformly 
continuous in x and measurable in f; see 1261 . 

If 1/ = R X Q, where the base Q is a bounded C' domain and if the coefficients A"^ of 
the operator ^ are functions in VMO v satisfying the conditions ( |1.2t and ( |L3t . then the 
condition (LH) is satisfied with parameters /io, A^i, and R,„ax depending on Q and UpiA^^) 
as well as on «, m, v. Therefore, we have the following corollary of Theorem |3.6l 
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(4.13) 



Corollary 4.11. Let 11 — ^ x CI, where Q is a bounded C' domain. Assume that the 
coefficients A°^ of ^ belong to VMO.v and satisfy the conditions ( 11.21 ) and (11.3b . Then, 
the Green 's matrix 0{X, Y) of I£ in 11 exists and for all X,Y e tl with X ^ Y and for all 
jiQ e (0, 1), we have 

mt, X, s,y)\ < Cx(0,ao)it - s) ■ 6iX, IT [(t - s) A RI,„\"'^ exp (-^^^^) > 

where d(X, Y) is defined as in ( 13. 7t . and R,„ax, C, and k are positive constants depending 
on Q and cjpiA"^) as well as on n, m, v, and hq. 

In the above corollary, one may assume that A"^ satisfy the weaker Legendre-Hadamard 
condition and may even include lower order terms in the operator. More precisely, let 

(4. 12) if^M = M, - Da{A"PDpu) + D„{B"u) + B"D„u +Cu + Au, 

where A"^, B", B", and C are m x m matrix valued functions on R"^' satisfying 

A"j'(X)^^^'r]p?]a>v\^f\Tjf, VfeR"', Vi/ e R" , VXeR"+'; 

f llA'^q' < y-'; f {\\B"\L + iml) + ||c||' < y-\ 

/ 1 M ML" / i \ll ML" M ML™ / M ML°° 

ff,/3=l a=l 

for some constant v e (0, 1], and A is a scalar constant. 

Corollary 4.14. Assume 'ZY = R x O, where Q. is a bounded C' domain. Let the operator 
Jff^ be as in ( 14.121 1 with coefficients satisfying the condition ( 14.131 ). We assume further that 
the leading coefficients A"^ belong to VMO.,. There exists Aq > such that if A > Aq, then 
the Green's matrix Q{X, Y) of S£x in tl exists and for all X,Y e tl with X ^ Y and for all 
jjLQ e (0, 1), we have 

mt, X, s,y)\ < Cx(0,oc)(t - s) ■ 6{X, Yf [{t - s) A RI,,\"'^ exp Lk^^^\ , 

where 6{X, Y) is defined as in ( 13.71 1. and Rmax, C, and k are positive constants depending 
on Q, u)p{A"^), and A as well as on n, m, v, and fiQ. 

We give a sketch of proof for Corollarv l4. 141 First we note that for sufficiently large A, 
one has the solvability of the following problem in the function space 'f^ ' ((■*, °°) x ^)'": 

u{s,-) = 0, 

where / e L'^i'U). In particular, one can construct the averaged Green's matrix ^{X, Y) 
of ^1 in 1/ by following the argument in ||8] §4]. We also note that the condition (LH) is 
satisfied in this case; see e.g., IfTTI . Then by modifying the proofs of Theorem |3.1| and |376l 
one can prove the above corollary. The details are left the the reader 

Remark 4.15. In Corollarv 14. 1 1 1 and Corollarv l4.14| the conditions of O and A"^ can be 
relaxed. We may assume that O is a bounded Lipschitz domain with a sufficiently small 
Lipschitz constant, and (jjp(A"^) is also sufficiently small for some p > 0; see e.g., IfTTI . 
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5. Proofs of Main Theorems 



5.1. Proof of Theorem 13111 By Lemma [Q and E Theorem 2.7], the condition (IH) 
implies the existence of the Green's matrix QiX, Y) of ^ in 1(. In fact, we point out that 
in the proof of fS] Theorem 2.7], one can completely replace the property (PH) by the 
condition (IH), the latter of which is weaker. This observation is useful because in the 
presence of lower order terms in the operator .if, the property (PH) does not follow from 
the condition (IH). Notice from fSl, Theorem 2.7] that we have 

Q{t,x,s,y) = if t < s. 

Therefore, to prove estimate ( I3.2l i. we only need to consider the case when t > s. To 
derive the estimate ( 13.21 ). we modify the proof in [S] §5.1], which was based on that in Il22il . 
We mention that the method in [22| was in turn based on the ideas appeared in ||9l and ifTSJI . 

Let i/' be a bounded Lipschitz function on R" satisfying \Dtfr\ < j a.e. for some y > to 
be chosen later. For r > i, we define an operator P ^^^ on L^(Q)"' as follows: For a given 
/ € L-(Q)"', let u be the unique weak solution in ^°'((i, oo) x Q)™ of the problem (see E 
Lemma 2.1]) 



f .iftt = 0, 
\ u(s,-) = e 



and define P'\.^J{x) :- e'^''-''^u{t, x). Then, by |[8!, Theorem 2.7], we find 

(5.2) PUf(x) ^ e^^-'> f g{t,x,s,y)e-^^y^f{y)dy, MfeL\Q.r. 

Jo. 

Then, as in |I8] §5.1], we derive 

(5-3) \\pU\U^ ^ -'^^^'-1|/|L.(a)' ^^ > ^' 

where i? = v"^ . We set p = y/t - s A R„iax and use the condition (LB) to estimate 

e-'^^-^yUf(x)f^\u(t,xf 

<A^,y"+2) r r \u(T,yfdydT 
Jt-p'- Jn,,(.v) 

where Slp{x) := D.n Bp(x). Thus, by using ( 15. 3t , we derive 

\PUf(x)\'<Nlp-"-' f f .2«-''-2^W|Pt,/(j)f^yJr 
Jt-p'- Jn,,(A) 

<N^,p-"-'- f f e^yP\pUj{yfdydT 

Jt-p- Jq,,(x) 

<Nlp-"-^e^yp f e^'^'^'-^^fWl^^^.^dr 
Ji-p- 

We have thus obtained the following L? — > L°° estimate for P^^,: 
(5-4) ll^^-/|L^(n)^^op-«^^.-^'^^<'-1|/|L.(„, 
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We also define the operator Q,^s ^^ L^ii^)"' for s < thy setting Qf^sSij) - e"'^*-'-' v(i, y), 
where v is the unique weak solution in ^ ' ((-oo, t) x D)'" of the backward problem 

'Jf V = 0, 
v{t, ■) = e'^g. 



(5.5) ^ ..,. , _ . 



By a similar calculation that leads to ( 15 .4t . we obtain 

(5.6) \\QUMl^,a,<Nop-"''e^^-^y-^'-%l,^^ 

Notice that it follows from (lOT l and (l53T l that 



fiPt.f)g- ffiQUg). 
Jii Jn 



Jii Jn 

Therefore, by duality, ( 15.6b implies that for all f,g& L^°°(Q)"', we have 

(5.7) lkt/IL,a)^M'P""^'-'"^^'"'"1|/IL.(n)- 
Now, set r = (s + t)/2 and observe that by the uniqueness, we have 

pLJ = PUiPUn V/ e C(Q)'". 

Then, by noting that t - r - r — s - (t — s)/2 and p/ Vz < '\/t - r A /?niax ^ p, we obtain 
from (Isil l and (|5J] | that 

ll^^r/|L^,a) ^ Cp-« ^^^'^^^^^^'-i/IL.a,, V/ 6 L^inr, 

where C - 2"^^A^y. For all x,y e Q. with x + y, the above estimate combined with ( I5.2l l 
yields, by duality, that 

(5.8) e*(-*'>-'^W|^(f^ ^^ j^^;)! < Cp-n g2yp+»yH,-s)_ 

Let (/'(z) := yi^odz ~" yl)- where (/tq is defined on [0, oo) by 

r if r < |jc — yl 



^o('-) , 

|x — yl if r > |x - yl. 

Then, i/r is a bounded Lipschitz function satisfying \D4r\ < y a.e. Take y = |x - 3'|/2)?(f - s) 
and set ^ := |:c - yl/ Vf- i. By ( 15.8b and the obvious inequality p/ y/t - s < 1, we have 

\0(t,x, s,y)\ < Cp-" exp{^/& - f/4^}. 

Let A^ = N(§) - N(v) be chosen so that 

exp(^/i? - f/4-id-) <Nexp{-f/M), V^ € [0, oo). 

If we set K - 1/8)? = v-'/8, then we obtain 

\git, X, s,y)\ < Cp-" exp [-k\x - y\-/(t - s)} . 

where C = C(«, m, v, No) > 0. We have thus proved the estimate (13.2b . ■ 



GREEN'S MATRIX FOR PARABOLIC SYSTEMS 



5.2. Proof of Theorein l3.3l As mentioned in the proof of Theorem FS . 1 1 the condition (IH) 
impUes the existence of the Green's matrix 0{X, Y) of ^ in TY. Moreover, the Green's 
matrix '0(X, Y) of 'Jff in 1/ also exists and we have the following identity: 

(5.9) g{X, Y) = 'g(Y Xf, \/X,Yet{, X+ Y. 

With aid of the above observation, we shall prove below that ( 13.4b implies part ii) in the 
condition (LB). Thanks to ( 15. 9b , the proof that ( 13.41 ) implies part i) in (LB) is similar and 
shall be omitted. Notice that ( 13.4b implies, via straightforward computation, that 

(5.10) \Q{X, Y)\ < C\X - Y\;;, if < If - .1 < Ri,,, 

where C = C{n, m, v, Co, k). Then, by the energy inequality (see [8] Eq. (3.21)]) and ( 12.21 ). 
and observing the obvious fact that Rma.J2 and R^a.^ are comparable to each other in the 
case when /?,„„, < oo, we obtain 

(5.11) \\g(; i')IL2.4/„(^\a.(y)) + \mi; miii^Q^^Y) ^ Cr-"l\ Vr € (0,i?„,,), 

where C = C(n, m, v, Co, k). 

Let X e 11 and R e (0, Rmax) be given. Without loss of generality, we may assume 
X -Q and for simplicity of notation, we shall write IJ^ - 1J^{Qi), etc. 

Assume that « is a weak solution of '^u = in 1/^ vanishing on S^. Let w - ^u, 
where ^ is a smooth cut-off function on R"^' satisfying 

(5.12) < ^ < 1, supp^ c Qri2, ^ s 1 on fofl/s, m\ < l^/R, and |f,| < 16/R^. 
Then w is a weak solution of 

'^w = -(,u - 'A"PD„(DpU - Da('A"PDpi;u) 

in 11^ :- R+ x Q.. Notice that w vanishes on R+ x 5Q and w = in {R~IA, oo) x Q.. 

For Y e K^/^, let p > be such that Qp(Y) c 1/+. For A: = 1, . . . , »i, let Vp be the A:-th 
column of ^(s Y), where 0^(-, Y) is the averaged Green's matrix of ^ in 'ZY as constructed 
in E §4.1]. Then, similar to lH Eq. (3.8)], we have 

(5.13) f ^"' = - r i,uvp- { (!A"PD,(Dpu)-Vp+ { (A^^Dp^u) ■ D,Vp 

-. Ii+h + h. 



which simply means that 



f W^Vp= f 'J^WVp 



Notice from ( |5l2] i that \Y - Z\p > R/d. for Z € 1/+ n suppD^, and recall 'A"^ = {A'^''f. 
Thus, if we set r - R/E A {dy A R^), then after interchanging indices a andyS, we find that 



Io + 



h^-f A';fg'{;Y)D,u'Dp(+ f A';fDpg'{;Y)u'D„^. 



Then, in light of HI Eq. (4.15)] and Id Eq. (4.16)], we take Umits p to zero in i5A3[ 



and conclude that for a.e. Y 6 'H^/4, we have 



(5.14) u'iY)^- f C,uJgj,{;Y)- r ^J^gjk{;Y)DjDp^ 

f A';^Dpgjk{;Y)u'D,C^:I[+I'^ + I'^. 



J-^m 
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Let JIr{Y) := 'U3r/4(Y) \ Qr/s(Y) d 'ZZ+^j \ Qr/s{Y). By (EH, Holder's inequality, and 
( 15. lit , we estimate 

|/;| < CR--\\g{; i')||i2(^^(j,)) ||H|b(«J,^) < CR-^"^^^'%h2^V.y 

Similarly, by ( 15.121 ) and ( 15.11b . we obtain 

I/3I ^ CR-'\\Dg(;Y)\y^^,^^r)) MkHvy) < CR-^"^'-^'^\\u\y^^.y 
Notice that the energy inequality (see e.g., 1281 §111.2]) yields 

(5.15) I|£>«IIl^(«;„) < CR-'\\uhny.y 
By ( I5.I2I 1. ( 15.1 Oi l, and (15.151 1, we estimate 

|4| < CR-'\\g{; i')|b(^^,,(y)) \\Du\y^t,;^^j < CR-^"^^^'^\\u\y^^,.^y 
By combining above estimates for I[,I'2, and /j, we conclude from (15.14b that 

where C = C(n, »j, v, Co, /c). Since the above estimate holds for all X € 1/ and R € (0, R„,ax), 
we obtain (LB) by a standard covering argument. ■ 

5.3. Proof of Theorem l3.6l Notice that by Lemma [Q] and Theorem l3.1l we have 

(5.16) \git,x,s,y)\ < CoXioMt-s) ■ [{t - s) ARi^,y"'\xp[-K\x-y\^/{t - s)] , 
where Co = C()(n,m,v,jiQ,Ni). We denote 

^i(X, F) = 1 1 A 1 and SiiX, F) = ( 1 A ^ 



Rmax f^\X -y\pl \ Rmax f^\X -y\ 

SO that 6{X, Y) - 6i{X, Y) ■ 62{X, Y). To prove the estimate dlSl ). we first claim that 
(5.17) \g{t, X, s,y)\ < C;r(o,oo)(f - s) ■ 6,{X, 7)"° {(r - .) A R;,,,,^"'^ exp {-^^3^ 

where C = C{n, m, v,jj.o,Ni). The following lemma is a key to prove the above claim. 

Lemma 5.18. Let 1/ = R x Q and assume the condition (LH). For R e (0,R,„ax) <^nd 
X & 14 such that dx < R/2, let u be a weak solution of J^u — in ll^iX) vanishing on 
Sg{X). Then, we have 

(5.19) \u{X)\ < C4°i?-"/2-'-'"'||H|b(^^-(x)), 
where C — C{n, m, v, fiQ, Ni). 

Proof. By the very definition the condition (LH), we have 

(5.20) \u(X') - u{X)\ < C\X' - X|7 ^-"/'-'-''°||h|Il2(«„-(z)), VX' € Q-r,,(X). 

For r e (dx, R/2), there is X' e g^/j W \ '^ ^^^^ ^^^ 1-^ - X'\p ^ r. By (l5^20l i we obtain 

|h(X)| = |h(X) - u{X')\ < Cr^'^R-'^'^-'-^'^My^^-^^xyy 
By taking limit r ^ dx in the above inequality, we derive (15.19b . ■ 

Now we are ready to prove ( 15.171 ). Take R - {R,„ax /\\X — F|p)/4. We may assume that 
dx < R/2 and t > s because otherwise ( 15.17b follows from ( 15.161 ). We then set u to be the 
A:-th column of g(-, Y), for k - 1, . . . , m, in Lemma l5.18l to obtain 

(5.21) \g(x, Y)\ < cd>'^"R-"'^-'-'">\\g(; Y)\y^y-^^x)); R = (R,nax a\x- Y\p)/4. 

Next, we consider the following three possible cases. 
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Case 1: \x - y\ < ^Jt - s < R,„ax- In this case R - Vf - ■s/4 = 1^ - Y\plA and thus, we get 
from dl^B and (ITIOl l that 

\G(x,Y)\ < cd^R-"^'-'-'"'\\g{;Y)\y(v-^^x^) < Cd'^^^R-"-"", 

which immediately implies ( 15.17b in this case. 

Case 2: -\/t - s < |x - yl A Rma.t- In this case R - (\x - y\ A R„,ax)/^- We denote Z = (r, z) 
and claim that for all Z e l/j^CX), we have 

(5.22) Mr, z, s, y)\ < CCo(f - s)-"'^ exp [-k\x - y\^/4{t - s)) , 

where Co and k are the same constants as in ( 15.161 ) and C - C{n, k). To prove the claim, 
first note that we may assume 7 = without loss of generality. Then by ( 15.161 ) we have 

mr,z, s,y)\ < CoX(0,^)(r) ■ r-"l^e-'^^'l'- < CoX(0,^){r) ■ r-"/2e-l-F/4.^ 

where we used |z| - \z - y\ > \x - y\l2 - |x|/2. Let us denote 

Then the claim ( 15.221 ) will follow if we show that there exists a positive number C - C{n, k) 
such that g{r) < Cg{t) for all r < f < \x[^, which in turn will follow if we show that 
go{r\) < Cgo(r2) for all r^ < r2 < I. But the latter assertion is easy to verify by an 
elementary analysis of the function go. 

We have thus proved ( I5.22l i. which combined with ( 15.21b yields 

\g(X, Y)\ < Cd''^°R-'">(t - s)-"'^ exp [-k\x - yf/4(t - s)] . 

Therefore, we also obtain ( 15.171 ) in this case. 

Case 3: Rmax ^ Vf - s. In this case R - Rmaxl^, and the desired estimate ( 15.171 ) becomes 

(5.23) \&{t,x,s,y)\ < C{dxlR,nnxrRA.^M'K^\x-y\^l^{t-s)]. 
Since t- s> 16R^, for all Z = (r, z) e V:^^{X), we have 

(5.24) exp \ -K \ < exp < -k \ < e ' exp ' 



r - s ) \ t — s J (^ 2(f — s) 

Then, from (15.21b . (15.16b . and ( 15.24b . we obtain ( 15.23b . which implies (15.17b in this case. 

We have thus proved that the estimate ( 15.17b holds in all possible cases. Finally, notice 
that the condition (LH) is symmetric between Jff and 'J^. Therefore, by repeating the above 
argument to 'Q{-,X), using the identity ( 15. 9b . and utilizing the estimate ( 15.171 ) instead of 
( 15.16b . we obtain (13.8b with a-/ 16 in place of a-. The theorem is proved. ■ 

6. Green's matrix for elliptic systems in two dimensional domains 
In this section, we are concerned with the Green's matrix for elliptic systems 

(6.1) Lu := -D„{A°l^{x)Dpu) 

in a domain Q c R^. Here, A"^ - A"^{x) are m x m matrix valued functions on 'E? with 
entries A. (x) satisfying the strong ellipticity condition 

(6.2) A",^{x)^l^>^>v\4\\ V^eR2'«, Vx g R^^ 
and also the uniform boundedness condition 



(6.3) ZZK^I'^^"'' ^-^^- 
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for some constant v e (0, 1]. We emphasize that we do not impose any other conditions on 
the coefficients. The adjoint operator 'L is defined by 

'Lu = -D„(A"I^Dpu), where 'A"P = {AP"f. 

Note that the coefficients 'A' of 'L satisfy the conditions ( 16.2b . ( 16.31 ) with the same v. 
Throughout in this section, we shall always mean by ^ the operator 5, - L on R^; i.e., 

^u^u,- Dc,{A"P{x)Dpu) 

for functions u defined on the cylinder 1/ = R x Q. It is obvious that the operator ^ 
satisfies the conditions ( 11.2b and ( 11.3b . Moreover, by ll25l Theorem 3.3], the condition (IH) 
is satisfied by ^ = (5, - L; see also [8l Corollary 2.9]. 

6.1. notation and definitions. For p > I and k a nonnegative integer, we denote by 
W*''^(Q) the usual Sobolev space; see e.g., Il20l . The function space Y^ (Q) is defined as 
the set of all weakly differentiable functions on Q. such that Du € L~{Q) and rju e W^ (Q) 
for any 77 e C^(R-). An open set Q c R- is said to be a Green domain if 

[Xnii-ue cr(R')) € wl;\si). 

We recall that if f2 c R^ be a Green domain, then Y^' (fi) is a Hilbert space when endowed 
with the inner product 



(m, v) ;= I DjuDjV. 
Jn 



Moreover, C^{Q.) is a dense subset in this Hilbert space; see e.g., 1321 §1.3.4]. 

For a given function / - (/',... ,f'")^ e L^^^iQ.)"', we call u = (m', . . . , m'")^ a weak 
solution in Fg'^CQ) of Lh = / if « e Y^--{Q.) and 

(6.4) f A'^fD/^uJD.^cp' ^ f/V, V0eC(Qr. 

Jn Jn 

It is routine to check that if fi is a Green domain and « is a weak solution in Y^^'''{Q.) of 
Lu - 0, then u = 0. Therefore, a weak solution in Yq' (Q) of Lh - f is unique. 

For a Green domain Q, we say that an m x m matrix valued function G{x, y), with entries 
G,y(x, y) defined on the set {(x, y) € OxQ : x ^ y), is a Green's matrix of Lin O if it satisfies 
the following properties: 

i) G{-,y) e w!-\iJ) and LG(-,y) = 6yl for all y € O, in the sense that 



I 



,ff/S 



-'n 

^1,2/ 



A';fDf,Gjt(;y)D„(/>' = ,^''(3;), V0 € C^iQ)"'. 



'J 



ii) G(-,y) 6 Fi'2(Q \ 5,(3;)) for all y 6 Q and r > 0, and G(-,y) vanishes on dO.. 
iii) For any / = (/',..., /'")^ 6 L^(Q), the function u given by 



hW:= f G(y,x)f{y)dy 
Jn 



belongs to F^' (Q) and satisfies 'Lh = / in the sense of ( 16.4b . 

By the remark made above, part iii) of the above definition gives the uniqueness of a 
Green's matrix if Q is a Green domain. We shall hereafter say that G(x,y) is "the" Green's 
matrix of L in a Green domain Q if G(x, y) satisfies all the above properties. We define 

£)(n):= Vi^A^(Q), 
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where \Q.\ the Lebesgue measure of Q and 

^(Q) :- inf{dist(/i, ^2) : ^ hes between two parallel lines Zi, /i). 
We remark that Q is a Green domain if g{Q.) < oo; see e.g., IfTS], Lemma 3.1]. 

6.2. Main result. Recall that we assume that the operator L in ( 16. Il l satisfies the conditions 
( I6.2I 1 and ( 16.3b . and that Q. is an open connected set in R-. In the sequel, we denote d_^ - 
dist(x, do.) and ln+ t - max(ln f, 0). 

Theorem 6.5. Assume that g = p(Q) < 00. 

(a) Suppose the condition (LB) is satisfied by ^ — d, — LinlJ — 'RxQ. Then the Green 's 
matrix G(x,y) of L in Q exists and we have 

(6.6) \Gix,y)\<cl—^]\cinJi^^], \fx,yen, x+y, 

XQ/^Rmaxj \ \x-y\ j 

where C — C{m, v. No). 

(b) Suppose the condition (LH) is satisfied by ^ — dt — L in tl — 'R X D.. Let G{x, y) be 
the Green's matrix of L in Q. Then for x,y &Q. with x + y,we have 

(6.7) \G{x,y)\ < c(l A ^r (1 A ^r [\ + In, ^^ ^ ^" 



V-y\) \ V-y\) \ \\x-y\ 



\QhRmaxl \ Q^RmaxI \ Q ^ Rn. 

where C - C(m, v,po,Ni). 

Remark 6.8. We point out that the assumption g = g(0.) < 00 in Theorem l6.5l can be some- 
how relaxed. In fact, the quantity g is related to the constant K - K{Q.) in the following 
Poincare's inequality (see lfT2l Lemma 3.1]): 

II^IIl^cq) < ^IIV^IL^(n), V^ e C(Q). 
We can replace ghy K in Theorem l6.5l as long as such a constant K < 00 exists. 

Theorem 6.9. Assutne that the condition (LH) is satisfied byJif — dt — Lintl — RxQ. 
with R,nax — °°- Then the Green's matrix G{x,y) of L in Q exists and for all x,y E Q. with 
X + y,we have 

( d, Y" I d, Y" ( Id, A d, 

(6.10) |G(jc,y)| < C 1 A --^ lA--^ \\+\n ' ' 



\x-y\l \ \x-y\l [ \\x-y\ 

where C - C(m,v,fj.o,Ni). 
6.3. AppUcations. Here are some easy consequences of the main theorems in this section. 

Corollary 6.11. Assume m — \ and g — g{^) < 00. Then the Green's function G(x,y) ofL 
in Q exists and for x,y &Q. with x +y,we have 

G{x.y)<C{\+\n^{gl\x-yW 

where C — C(v). If in addition, Q satisfies the condition (S), then for x + y,we have 

G(.,,).cfiA^rfiA^r(i-in.^^^'" 



\x-y\) \ \x-y\) \ \\x-y\ 



.c\^JU.J^VU. "' "■ 



[g A Ra j \ g A Ra j \ g A Ra 
where C - C(v, 0). 
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Proof. Notice that the condition (LB) is satisfied by ^ = 5, - L in any cyHnder 1/ = R x Q 
with Rmax - °° when m - 1; see e.g., Il30l Theorem 6.30]. Moreover, the condition 
(LH) is satisfied as well with R,„ax - Ra if the base Q enjoys the condition (S); see Il30l 
Theorem 6.32] and Lemma [8^ Therefore, the corollary follows from Theorem l6.5l ■ 

Corollary 6.12. Let Q. — [x &'Rr : X2 > ip{x\)], where (^ : R ^ R « a Lipschitz function 
with a Lipschitz constant K. Then, the Green's matrix G{x,y) ofL in Q exists and for all 
x,y & Q. with x +y,we have 

I d, r / d, r f (d, A d. 



\x-y\j \ \x-y\j \ \\x-y\ 

where C — Ciin, v, K) and hq — jJLoiy, K). 

Proof. It is known that the condition (LH) is satisfied by^ = 5, -LinTY^RxQ with 
Rmax - °°; see IfTSl Lemma 4.4]. Therefore, the corollary follows from Theorem |6.9l ■ 

Corollary 6.13. Assume that Q is a bounded Lipschitz domain. Then the Green 's matrix 
G{x,y) of L in Q exists and for all x,y 6 Q with x +y,we have 

|G(x,y)|<c(lA^^) [\K-^lJ\ (l+ln,f ^ 



\x-y\) \ \x-y\) \ \\x-y\ 

where C — C(m, v, Q). 

Proof. By using |12] Lemma 4.4] and standard partition of unity argument, one can show 
that the condition (LH) is satisfied hy .5f - dt - L. The corollary is an easy consequence 
of Theorem l6.5l part (b) together with the assumption diam O < oo. ■ 



6.4. Proof of Theorem 16.51 Existence of the Green's matrix G{x,y) of L in Q follows 
from the assumption p(Q) < qq; see llT2l Theorem 2.121. To be more precise, first we recall 
that (IH) is satisfied by ^; see f25\. Let ^(f, x, i, y) be the Green's matrix of ^ in 1/ given 
as in Theorem l3.1l Since A"^ are independent of f, we have 0{t, x, s,y) - Q(t - s, x, Q,y). 
In the sequel, we shall denote the "Dirichlet heat kernel" of L by 

K(t,x,y) = K{X,y) := g{X,Y) = g(t,x,Q,yy, Y = (O,^). 

It is shown in |[T2l that if £) < oo, then the Green's matrix G{x,y) of L in Q is given by 

(6.14) G{x,y)= I K{t,x,y)dt, 'ix,yeQ., x+y. 

Jo 

Part (a). From the estimates ( 15.10b and ( IS.lll l, we obtain 

(6.15) \K{X,y)\<C\X-Y\-/ if < t < rI,/, 

(6.16) \\K{;y)\\L2..i..^^^Q^^Y)) + lll^(-'3')lll«\e^(?) ^ '^'""'' ^'' < ''^"'- 

By using (LB) and (I6.16I) . and following the proof of lfT2l Lemma 3.12], we also obtain 

(6.17) \K{t,x,y)\<Cgr-^e-^'''^'-^''"''^\ Vf > 2r\ Vr e (0,R,„ax)- 

Now, we are ready to prove the estimate ( I6.6I 1. We set r = {gARmax)/'2.. If < \x-y\ < r, 
then by ( 16.141 ). we have 

(6.18) \G{x,y)\< \ +\ +1 \K(t,x,y)\dt-.h+h + h. 

Jo J\x-v\- Jlr- 
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It then follows from (16.151) and (16.171 ) that 

/i <C \x - yl'^ dt < C, 

Jo 

h<C \ r^dt<C + C \n(r/\x - y\), 

l3<C I Qr-\-^''^<-^''^le' dt < Cr-^g\ 
Combining the above estimates together we obtain that 

|G(.,3')l<c(^^f+Cln^^^^"- 



\g ^Rmaxl \ 2\x-y\ 

which proves the estimate ( 16.61 ) in the case when < \x - y\ < (g A R„,ax)/2. 

On the other hand, if \x - y\ > r = (g A R„,ax)/2, then we estimate by ( 16.151 ) and ( 16.171) 

p2r^ poo 

\Gix,y)\< + \K(t,x,y)\dt 

Jo Jlr- 

n2r- poo 

<C r-- + C er-3e-2H'-2^^)/e^ dt<C + Cr-'g' . 

Jo Jlr- 

Therefore, we obtain the estimate ( 16.6b also in this case. 

Part (b). By Lemma [STI and Part (a) above, we obtain the existence of the Green's matrix 
G{x,y) of L in Q and the estimate ( I6.6I 1 with C - C{m,v,iio,N]). By Theorem 13.61 and 
observing that 2R„,ax and /?,„„, are comparable to each other when R,„ax < '^, we find that 
if < |X - Y\p < 2Rmax, then we have 

/ d Y" I d V" 

(6.19) |»r(X,,)|sc(.A^) (lAj^) IX -?|- 

We claim the following estimate holds: For all r e (0, R,„ax), we have 

(6.20) \K(t, x,y)\ < Cgr^l A (dJr)Y"'{l A (dy/r)}'"'e-^''^'-*''^'<^\ Vf > 4rl 

Let us assume the claim (I6.20l i for the moment and prove the estimate (16.7b . Similar to 
( 16.18b . in the case when < \x - y\ < r :- {g A R,„ax)/2, we estimate 

|G(x,y)|< ' + + \K{t,x,y)\dt^:h+h + h, 

Jo J\x-y\- J4r- 

It follows from ( 16.191 ) that 

/ d, r/ d, ^^° 

/i<C lA--^ lA 



\x-y\l \ \x-y\ 



Po I A \fo 



/2<CllA- — ^1 |1A--^1 -Jl+ln 



r 



\x-y\l \ \x-y\l [ \\x-y\) 

Also, by ( 16.20b . we obtain 

/3 < Cg'r-'il A idx/r)r{l A (dy/r)}"". 
Combining the above estimates for /i,/2, and I3 together, we obtain ( 16.71 ) in this case. 
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lf\x-y\>r = {gA Rmax)/2, then by using (16.191) and (16.201) . we estimate 

p4/'^ poo 

|G(.^,3')I< + \K(t,x,y)\dt 

Jo J4/-2 

< C{1 A (djr)r{l A (d,/r)r + Cg'r'll A (djr)r{l A (^v/r))"" 
<Cg\-'{lA(dJr)r{lA(dy/r)r. 

Therefore, we also obtain ( 16.71 ) when \x - y\ > [g A R„,ax)/2. 

It only remains for us to prove (16.20b . The strategy is similar to the proof of ( 13.81 ). We 
first prove the following estimate, which is a "half" of the estimate (16.20b . 

(6.21) \K(t,x,y)\<Cgr-^lA(dJr)re--'^'-''-'^/e\ \/t > 3r\ Vr e (0,/?„,„). 

To prove ( 16.21b . it is enough to assume that d^- < r/2 because otherwise ( 16.21b follows 
from (16.17b . By Lemma 13 . 1 8 1 applied to k-th column of K{-,y), for k = 1 , . . . , m, we have 

\K(t,x,y)\ < Cfiff r-'-'^«||/ir(-,y)|b(i/,(x)), Vf > Sr^, Vr 6 (0,7?,„„). 

Then by ( 16.171 ). we obtain 

\K{t, x,y)\ < Cgd'^^r^-^'^e-^^^'-^''^!'^, Vf > Sr", Mr 6 (O,7?™,0- 

Therefore, we proved the estimate ( 16.211 ). Next, let us recall from lfT2l that the Dirichlet 
heat kernel 'K{t, x, y) of the adjoint operator 'L is given by 

'K{t,x,y)^K{t,y,x)^. 

By repeating the above argument to 'K{-,x), using the above identity, and utilizing the 
estimate (16.211) instead of ( 16.17b . we obtain ( 16.20b . This completes the proof of Part (b). 

The theorem is proved. ■ 

6.5. Proof of Theorem I6.9[ We use the same notation used in the proof of Theorem l6.5l 
Since Rmax - °°, we find that ( 16. 19b valid for all f > and x eQ.. For x,y e Q. with x + y, 
let r > max(|jc - y\, d^ A dy) be a number to be fixed later We write 

pea /^l-t-Vp pr- pea 

\K(t,x,y)\dt^ + I + I \K{t,x,y)\dt-.h+h + h. 

Jo Jo J|a-v|- Jr^ 

Then it follows from ( 16. 19b that 

h <C|1 A— ^1 (l A ■ 



l-'^ -y\l \ \x- y\ 

d Y° I d Y" 

/2 < C 1 1 A — -^- 1 A — -— ■ 2 In 



\x-y\j \ \x-y\j \\x-y\ 

ri-'^°^2 ^f ^ 2C(d^ A dyfiu^^W-f", 

where C = C{m, v,//o, M). 

We define G{x,y) by the formula ( 16.14b . It is shown in the proof of fT2|, Theorem 2.21] 
that thus defined function G{x,y) is indeed the Green's matrix of L in Q. 

Now, we will choose r as follows. 

Case 1: \x - y\ < d^ A dy. We take r - d^ A dy. Then 

\G{x, y)\<Ii+h+h<C + C\n [d,, A dy/\x - y\) + C, 
and thus, (16.10b is satisfied in this case. 
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Case 2: dx A dy <\x -y\ < dx V dy. We take r -\x - y\. Then we have 

/ d Y° I d '^^° 

|G(x,3;)|</i+/3<C lA^^ lA -' 



\x-y\i \ \x-y\, 

It is clear that ( 16.101 ) is satisfied also in this case. 

Case 3: dx V dy < \x - y\. We take r - \x - y\, and obtain a better estimate for I^. 






Therefore, we have 



d Y" I d ^''° 

|G(jc,y)|</i +/3 <Cll A--i— lA ■ 



\x-y\j \ \x-y\j 
and thus ( 16.101 ) follows in this case too. 

The proof is complete. ■ 

7. Global estimates for systems with Holder continuous coefficients 

In this section, we assume that the coefficients A"j in ( II. lb belong to the class ^i^/^-f, 
where ju € (0, 1 ), and obtain improved global estimates for the Green's matrix in '1/ = RxQ, 
where fi is a bounded C''' domain. We remark that these assumptions can be relaxed; see 
Remark lTTSl Recall that we have the following gradient estimate; see e.g, |1T9] Chapter 3]. 

Lemma 7.1. Under the assumptions above, there exist R,„ax 6 (0, oo] and C > so that 
for all X e tl — 'Rx D. and < R < Rmax, the following is true: Ifu is a weak solution of 
Jffu — in t(^{X) vanishing on »S^(X), then we have 

(7.2) I£>«l0;t/,,,(X) < CR-^-"'^\\u\\l2^^,-^^x))- 

Here, the constants R,„ax and C depend on n, m, v, A"^, and Q.. 

The following lemma is an immediate consequence of Lemma l7.1l c.f. Lemma l5.18l 



Lemma 7.3. Assume the conditions ofLemma \7.1\ For R e {0,Rmax) andX e 11 such that 
dx < R/2, let u be a weak solution of J^u — in ll^iX) vanishing on S^iX). Then, we 
have 

\u{X)\ < CdxR-^-"^^\\uh2m,(x)), 
where the constant C depends on n, m, v, A"^, and Q.. 
Now we state the main result of this section. 

Theorem 7.4. Assume the conditions of Lemma \7.1\ Let Q{X, Y) — Q{t, x, s, y) be the 
Green 's matrix of J^ in 11 = R X Q and let 6(X, Y) be defined as in ( I3.7l l. Then for all 
X,Y eK with X it Y, we have 

(7.5) IQit, X, s,y)\ < Cx(0,oc){t - s) ■ 6iX, Y) [{t - s) A Ri.xY"'^ exp (-^^^) , 

where k — k{v) and C is a constant depending on n, m, v, A"^ , and Q. 



Proof. We note that Lemma 17.11 implies the condition (LH) with jiq - 1 , which in turn 
gives the existence of the Green's matrix 0(X, Y) of ^ in K. Now we can repeat the 
argument in the proof of Theorem |3.6| bv using Lemma 173] instead of Lemma |5.18l ■ 
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Remark 7.6. Following 11291 Section 5], the same gradient estimates ( 17.21 ) is available for 
systems with Dini continuous coefficients in Dini domains, which is more general than the 
Holder conditions; see Il27ll for the definition of Dini domains. Since we only used the 
gradient estimate ( 17.21 ) for weak solutions in the proof of Theorem 17 .41 the conclusion of 
the theorem is still available for systems with Dini continuous coefficients in Dini domains. 
We would like to hereby thank YanYan Li for helpful discussion regarding the gradient 
estimates for systems with Dini continuous coefficients. We remark that a weaker form of 
the estimate ( I7.5l l for scalar Green's functions was obtained in Q by using the maximum 
principle argument. 

8. Appendix 

Lemma 8.1. Let 1/ = R X Q and assume the condition (LH). Then, the condition (LB) is 
satisfied with the same R,„ax ond No — No(n, m, /do, Ni ). 

Proof. We shall only prove i) since the proof of //j is very similar Let « be a weak solution 
of ^u = in tl^iX) vanishing on S^{X), where X elA and R e {Q,R,„ax)- By using the 
triangle inequality, for all Y e Qj^niX) and Z e 2r/2(^)' we have 

\u{Yt < 2\u{Y) - u{Zf + 2\u(Z)f < CR^'"'[a]l/2,,,;Q-^fx) + C|fi(Z)p. 
Then by taking average over Z e Q^niY) and using (LH), we obtain 



u 






where C - C(n, m,//o,A^i). The proof is complete. ■ 

Lemma 8.2. Assume that the operator J^ is given as in ( 12. Il l and satisfies the conditions 
il.li and ( 11.31 ). Then the condition (IH) is equivalent to saying that the operator J^ and 
its adjoint LSf satisfy the property (PH) in |[8l. 

Proof. We shall only prove that part i) of (IH) is equivalent to the property (PH) for ^ 
since the proof that part ii) of (IH) is equivalent to the property (PH) for '.if is very similar. 
Let H be a weak solution of ^h = in Q^(X), where X e 1( and < R < R^ A dx- 
Recall that the property (PH) in H] for ^ is as follows: We have 

(8.3) r \Duf- <Co(-T''"' f \Du\", W<p<r<R. 

jQ-piX) ^rl Jq-(X) 

First, we assume that the condition (IH) holds and prove ( 18. 3t . We may assume that 
p < r/4; otherwise, ( 18. 3t is trivial. We denote 

(")e^ = T "■ 

JQ7(X) 

Notice that we may assume, by replacing « by « - (u)q-, if necessary, that (h)q- - 0. 
From the energy inequality (see e.g., 1281 §111.2]), the condition (IH), and then parabolic 
Poincare's inequality (see e.g., ||8] Lemma 2.4]), it follows (recall thatp < r/4) 

r \Du\-<Cp--f \u - {u)q- \" < Cp-- f -f \u(Y)-u{ZfdZdY 

< cp-^-[«]2„,,,„^ „ < cp-^-.-^- £ |„|2 

•J Qi (X) 



Kcr-T'" I idhP 



JQ-(X) 
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where C = C(n, m, v,/io, Co). We have derived the property (PH) for J^. 

Next, we assume the property (PH) for J^ and prove the condition (IH). By the para- 
bolic Poincare's inequality (see e.g., [[8] Lemma 2.4]) and the property (PH), and then the 
energy inequality, we obtain for all Y e Q^.^{X) and r e {0,R/4-] that 



r \u-(u)Q-f<Cr' f |DH|2<Cr2(^)"^"'° f \Duf 
Jq7(Y) Jq;(Y) ^'<' JQr/AY) 

where C = C{n, m, v,^q, Co). Then by the Campanato's characterization of Holder contin- 
uous functions (see e.g., ||8] Lemma 2.5]), we obtain 

Then, the above inequality together with a standard covering argument yields part i) of the 
condition (IH). The proof is complete. ■ 



Lemma 8.4. Let t( = M.X Q and assume the condition (LH') in Remark \3.9\ Then the 
condition (LH) is satisfied. 

Proof. We shall only demonstrate that part i) of (LH') implies part i) of (LH) since the 
proof that part ii) of (LH') implies part ii) of (LH) is very similar. 

Let H be a weak solution of J^u = in tl^iX) vanishing on S^{X), where X etl and 
R 6 {<d,R„u,x)- Notice that from (LH') we have for all Y e QT^i^iX) and r e (0,7?/4] that 

(8.5) r \Du\^<c(!-y r \Du\\ 

where C = 4"^^''°Co. By Lemma |876l below, the above estimate ( 18. 5b , and the energy 
inequaUty, we have for all Y e Q^,^(X) and r e (0,R/4] that 

r \u- (u)q-\- < Cr^ f \Du\^ < Cr^ (^)"^"'° f \Du\^ 
Jq;{Y) Ji^;{Y) '^'<' JiJ-^iiiY) 

I f \n+2+2//o r r 

<C - \u\^ KCr^^^^^f^R-^"" f \u\-, 

where C = C(n, m, v,fio, Co). Then by the Campanato's characterization of Holder contin- 
uous functions (see e.g., H] Lemma 2.5]), we obtain 






Jq-ax) 



Then, the above inequality together with a standard covering argument yields part i) of the 
condition (LH). The proof is complete. ■ 

Lemma 8.6. Let 1/ = R x Q and u be a weak solution of ^u — f in 11^ iX) vanishing on 
S;{X), where f 6 L'"{V;{X)). Then we have 

(8.7) r |H-(H)e,f <Cr2 r |DH|2 + Cr2-"||/||2^ ; u ^ x^l,(X)U, 

•JQjiX) J'U^iX) 



where {u)q; - Jq-.v) " '^"'^ ^ - C(n, m, v). 
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Proof. We modify the proof of Ii37i Lemma 3] . Without loss of generality, we may assume 
X -Q. Let ^ = ^(x) be a smooth function defined on W such that 

< ^ < 1, suppf c B,., ^ s 1 on B,./2, and \D(\ < Air. 

Let 5"' = jg ^ > cr", where c = c{n) = 2""|Bi|. We denote 

Pit) -.^5 I ^{x)u{t,x)dx^6 I ({x)u{t,x)dx; j8 := r ^ | p{t)dt. 
By following the proof of l!37l Lemma 4], we estimate, for -R^ < s < t <Q, 



-2nii/'||2 

IIlI(«,t)- 



(8.8) \Pit)-p{s)\^ <Cr-" f \Du\- + Cr- 

Jtl; 

Since (h)^- minimizes the integral J„-\u - c\~ among c e R'", we obtain 

r |H-(H)e-i'< r i«-)8|'<2 r iM-^of +2 r ^sco-ysf. 

jQr Jq; Jq7 Jq; 

Notice that « e W,' (2,7) and Dm = xv; Du in 27- Therefore, by a variant of Poincare's 
inequality, we have 

f \u-p(t)\^^ f f \uit,x)-p{t)\^dxdt<C f \Du\-. 

Jq7 J-r- Jb, JllT 

On the other hand, by the definition offi and ( 18.8b . we obtain 

r - 9 r" - T r" r" 

|j8(0-j8| =|B,I |;8(0-)8r<Cr"-M |j8(0-j8(^)| flf^flff 

Jgr J-i- J-r- J-r- 

<Cr^ f \Duf + Cr--"\\f\\l^^^^. 
By combining the above inequalities, we obtain (18.71 ). The proof is complete. ■ 



Lemma 8.9. Assume a^^iX) satisfy the condition (14.71) . Let 1( -RxQ., where Q satisfies 
the condition (S), and define S as in ( 14. 81 ). where A^.{X) are the coefficients of the operator 
££. There exists So = (Soin, vq) > such that if S < Sq, then the condition (LH) is satisfied 
with fiQ - //o(n, vo, 0), Rmax — Ra, and N\ — N\{n, m, vq, 9). 

Proof. In this proof, we shall only consider part i) in the condition (LH) because proof of 
part ii) in the condition (LH) will be almost identical. 

We shall prove below that there is a positive number Si = Si{n, vo) such that if ^ < 
Si, then the following holds: There exist positive constants fii = jUi(n, vo,0) and Ci = 
Ci (n, m, Vo, 9) such that for any X € d^^ = R x dH and R € (0, Ra), if « is a weak solution 
of ^u = in tl^QC) vanishing on S^{X), then we have 

(8.10) r \Du\^ < Ci{-T ''' r \Du\-, ^0<p<r<R. 

We remark that the above condition is a paraboHc analogue of the property (BH) in ||231 . 

We also note that by [8 , Lemma 2.2], there is S2 - Szin, vq) > such that if S < S2, then 
the property (PH) in 18] is satisfied by ^ with the exponent 1J.2 = fJ^iin, vq) and Re = 00. 
More precisely, if « is a weak solution of ^u = in Q]f{X) c 1/, then we have 

(8.11) r |Dh|2 < C2 (-]"'" I \Du\~, VQ<p<r<R, 
jQp(X) ^rl Jq-^x) 
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where C2 = C2(n, m, vq). Then we combine (18. 10b and ( 18. Ill ), via a standard method 
in boundary regularity theory (see e.g., lfT9l §3.4]) and Lemma 18741 to conclude that if 
(ff < (ai A S2 ='■ ^a, then the condition (LH) is satisfied with parameters juq - fj.i A H2, 
N\ - Ni{n, m, vq, 6), and R„,a.^ = Ra', see Remark [33] 

By the above observation, it only remains for us to prove the estimate (l8.1Qb . For 
X e dp^ and R e (0, Ra) given, let « be a weak solution of ^u = in 1J^{X) vanishing 
on »S^(X). Denote by ^q the parabolic operator acting on scalar functions v as follows: 

J^ov = V, - DJa"f'Dpv). 

For r e (0, R\, let v' be a unique weak solution in 1^('Z/,7(X)) of the problem 

S^ov' = in ^l-{X)■, 



v' = u' on dpti;{X), 

where ; = 1, . . . , m. We claim that there are positive constants fj. = fi(n, vq, 0) and C - 
C(n, m, Vq, 6) such that the following estimate holds: 

(8.12) r \Dv\^ <cl-X'' f \Dv\^, VO<p<r. 

Notice that we may assume thatp < r/8 because otherwise ( 18.121 ) becomes trivial. Since 
each v' vanishes on S^{X), it follows from ll30l Theorem 6.32] and 1301 Theorem 6.30] that 
there exist yu = //(«, vo, 6)> Q and C - C{n, vq, ff) > Q such that 

(8.13) osc v' < Cffri" sup Iv'l < Cffr 



-j"-n/2-l||,,i 



V 



Il^cm-^c*))- 



In particular, the estimate ( 18. 13b implies v'{X) - 0. Then, by the energy inequality and 
Lemma 4.2], we obtain (recall that p < r/8) 

r |Dvf< Cp-2 r \vf^Cp-^{ \\KY)-v\X)fdY 

JnApiX) •J'U^AX) J'U^^iX) 



\cr-r'r--l Ivf 



iT'-i,,- 



< Cp" osc v' 

\'U^_„(X) 

<C(^| \Dv'\-, i^l. 



'^ri 



where C - C(n, vq, 0). This completes the proof of the estimate ( 18.12b . 

Next, notice that w :- u -v belongs to 'fzi'U^iX)), vanishes on dp^^{X), and satisfies 

Therefore, by the energy inequality we obtain 

(8.14) r \Dwf-<CS'^ f \Duf-, 

Jti;(X) J«-(X) 

where S' is defined as in ( 14.81 ). By combining ( 18.12b and ( 18.14b . we obtain 

I \Du\- < C (-X^''' f \Du\^ + C<g-\ \Du\^, VO<p 

Now, choose a number jUj e (O,/!). Then, by a well known iteration argument (see, e.g., 
llTSl §111.2]), we find that there exists (ai such that if d" < <?i, then we have the estimate 
( 18. 10b . The lemma is proved. ■ 



< r. 
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